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1 Current status of Material Design Using Correlated Electron 
Systems 

1.1 Weakly Correlated Systems 

The dream of accelerating the discovery of materials with useful properties using 
computation and theory is quite old, but actual implementations of this idea are 
recent flTJ [2) . Successes in material design using weakly correlated materials, are 
due, to a large degree, to a two important developments: 

a) approximate implementations of first principle methods, which are relatively ac- 
curate and computationally efficient 

b) robust implementation of algorithms which are highly reproducible and widely 
available in well tested codes. 

Density functional theory based approaches give reliable estimates of the total en- 
ergy, and are an excellent starting point for computing excited state properties 
of weakly correlated electron systems. These approaches allows the evaluation of 
transport coefficients using very limited, or no empirical information, and are be- 
ginning to be used in conjunction with data mining techniques and combinatorial 
searches. 
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2 Kristjan Haule and Gabriel Kotliar 

1.2 Strongly Correlated Electron Systems 



Since a large number of interesting physical phenomena, such as high temperature 
superconductivity and large Seebeck coefficients, are realized in strongly correlated 
electron systems, there is a great interest in the possibility of carrying out rational 
material design with correlated materials IDG]. 

The theoretical situation in this area, however, is a lot more uncertain. For ex- 
ample, the issue of whether the two dimensional one band Hubbard model supports 
superconductivity or not is still very open fl2l . Given that this model is an extraor- 
dinary oversimplification of realistic materials, it is hard to contemplate explaining, 
let alone predicting experimental results in materials that require a much more elab- 
orate models for their description. The prospect of predicting properties of materials 
which have not yet been synthesized is even more daunting. In this chapter we will 
argue that this assessment is overly pessimistic, and we will give some reasons why 
we expect a rapid progress in the coming years through the interplay of qualitative 
reasoning, new theoretical and computational methods, and experiments. We will 
then describe some attempts to gain experience in this field, and some lessons that 
we have learned using thermoelectric performance as an example. 

1.2.1 Dynamical Mean Field Theory 

The advent of Dynamical Mean Field Theory (DMFT) removed many difficulties 
of the traditional electronic structure methods. DMFT describes Mott insulators, 
as well as correlated metals. DMFT combines ideas of quantum chemistry, such 
performing a full configuration interaction calculation (at a local level to avoid size 
consistency problems), and physics, such as carrying out a diagrammatic expansion 
around the band limit. DMFT treats quasiparticle bands and Hubbard bands on the 
same footing, and, unlike simpler approaches such as LDA+U, is able to describe the 
multiplet structure of correlated solids. The latter is being inherited from open shell 
atoms and ions. DMFT has been successful in accounting for the behavior observed 
in correlated materials ranging from plutonium to vanadium oxides and have even 
made some predictions, which have been successfully confirmed by experiment|8, 
?]. This suggests that the approach is reasonably accurate, in the sense that it gives 
a zeroth order picture of correlated materials not too close to criticality. Ten years 
ago, a combination of DMFT with electronic structure methods LDA+DMFT, was 
proposed J5] [TT] [10 1 and accurate implementations are being actively developed 
across the world. Just like LDA, these tools connect the atomic positions with the 
physical observables using very little information from experiment, and therefore 
they have the potential to accelerate material discovery. 

Predicting the phase diagram of strongly correlated materials is an extremely 
difficult problem. Correlated materials have many competing phases, which are very 
close in energy. This poses serious difficulties to direct numerical studies of model 
Hamiltonians because omitting small terms which are present in the Hamiltonian 
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of the actual material or finite size effects connected to boundary conditions can 
exchange the stability of two very different phases. 

DMFT divides the solution of the many body problem of a solid state system into 
two separate and distinct steps: the study of the evolution of the mean field solutions 
as a function of parameters and the computation of total energies for each DMFT 
solution. Common to many mean field approaches, a given Hamiltonian can have 
many distinct DMFT solutions, describing various possible phases of a material. 
Which phase is realized for a given value of parameters (temperature, volume, stress, 
doping concentration of impurities, etc.) is determined by comparing the free energy 
of the different DMFT solutions. A lot of important information can be obtained 
from the first step alone, when combined with experimental information. If one 
knows that for some value of parameters certain phase is realized in material, one 
can use DMFT to explore the properties of that phase, and optimize desired physical 
property, sidestepping the difficult issue of the comparing the free energies of the 
different competing phases which can be done at a later stage. The rational material 
design process should then suggest modifications of the material to stabilize the 
phase with desirable properties. 



2 The process of rational material design 

Figure[T]describes schematically the rational material design process. It begins with 
a qualitative idea, which is then tested by a calculation. One of the major advances 
of realistic DMFT implementations such as LDA+DMFT or GW+DMFT is that 
now this calculation can be made material specific, resulting in a set of predictions 
that can be tested experimentally . The experimental results can either rule out the 
qualitative idea, in which case the process stops, or reinforce and refine the idea. 
Experiments also help calibrate the computational methods, which in turn lead to an 
improved material specific prediction in the next iteration. The expected result form 
this process are materials with improved properties M\, M2 , M3, •. In addition, this 
process tests theoretical ideas in an unbiased way, deepens our understanding of ma- 
terials physics, and refines the accuracy of computational tools. Large databases of 
existing materials are being created (for example http://icsd.ill.eu/icsd/index.html), 
which are starting to be used, in combination with the first principle methods, 
for data mining techniques [?]. Using the crystal structure information from the 
database, the first principles methods can identify potentially promising materials, 
which can then be analyzed experimentally. 
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3 Thermoelectricity of Correlated Materials 
3.1 Formalism 

The transport coefficients that govern the thermopower, electric and thermal con- 
ductivity can be expressed in terms of the matrix of kinetic coefficients A,„ relat- 
ing the electric and thermal currents J, Jq to the applied external fields V/i/T, 
VT /T 2 . Transport quantities become S = — (&;?/e)(Ai/Ao), <7 = (e 2 /T)Ao, K = 
k]j[A2 — Ai 2 /Aq\. The thermoelectric response thus reduces to the evaluation of ki- 
netic coefficients. 

The thermoelectric figure of merit is defined by 




Fig. 1 A schematic drawing of the rational material design process. It relies on condensed matter 
theory, material databases and realistic DMFT implementations, and it involves a close and iterative 
interplay of theory and experiment. 
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where T is the absolute temperature, a is the electrical conductivity, S is the Seebeck 
coefficient or thermopower, and K (K p i, (mon ) is the electron (phonon) contribution to 
the thermal conductivity. 

The Wiedemann Franz law is an approximate relation that allows us to estimate 
the ratio of the electronic contribution to the thermal conductivity (k) and electric 
conductivity (a). It postulate that the Lorentz number, L—k/ (oT), is weakly ma- 
terial dependent. 

Its value at low temperatures is given by (n 2 /3)(k B /e) 2 = 2.44 x 10~ S WQ /K 2 . 
We will return to the Lorentz number at higher temperatures later in this article. If 
we ignore the thermal conductivity of the lattice, the figure of merit can be written 
as ZT = S 2 /L, hence to have a promising figure of merit (ZT close to or larger than 
one) it is necessary to have S bigger than the basic scale k/e = 8610 6 V /K. The 
thermal current of an interacting electronic system was determined first by Mahan 
and Jonson [13 |. Ref ifTSIl discusses a model containing electrons interacting with 
phonons, and the review [18| discusses the general case of the electron electron 
interactions (see also ref [24Q. 

DMFT expresses the one particle Greens function in terms of a local self energy 
of an impurity model, satisfying a self consistency condition. Practical evaluation 
of the transport coefficients becomes possible in the approximation on small vertex 
corrections. This was first done by H. Schweitzer and G. Czycholl [27 1 (see also 
ref E5\ ). For the Hubbard like interactions, there are no contributions from the 
non local Coulomb interactions, and the negligence of the vertex corrections can be 
justified rigorously in the limit of infinite dimensions 1 15 1. The same is true, but far 
less obvious, for the thermal current, as it was shown in ref [24|. In the multiorbital 
situation, the vertex corrections to the conductivity need to be examined on a case 
by case basis, and do not necessarily vanish, even in infinite dimensions. With this 
approximation, the LDA+DMFT transport coefficients reduce to 

A£ v = itrfdco (^y^iv^p^Ki^ico)} (2) 

where Vk = — j| ( y/ki | V | y/k j ) are velocities of electrons and pt is the electron spec- 
tral density 

Pk(ffl)=2^[c£(ffl)-Gk(fl>)]. 0) 

The weakly interating case appears as a limiting case where the spectral function 
becomes a delta function pk,;(e) = 8(e — £ki)8jj. One can therefore formulate the 
problem of the optimization of the figure of merit as the problem of optimizing 
a functional of spectral functions, with self energies which are realizable from an 
Anderson impurity model, with a bath satisfying the DMFT self condition. 
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3.2 Thermoelecricity near the Mott transition: qualitative 
considerations 

Following the early developments of DMFT and its successful application to the 
theory of the Mott transition in three dimensional transition metal oxides (for re- 
views see [8 1 [?]), it was natural to use this approach to formulate and answer the 
question of whether we should look for good thermoelectrics near the Mott localiza- 
tion derealization transition. This talk describes our current understanding of this 
issue and the tentative answer, at this point, is no, but perhaps yes. 

There were several reasons to suspect that proximity to the localization dereal- 
ization transition is good for thermoelectricity: 




1 10 100 1000 

y (mJ/ mol K 2 ) 

Fig. 2 Bhenia Jaccard Flouquet plot from 1 4 1 . The theoretical point obtained on the LaSrTiCB 
system with twenty percent doping away from the Mott insulator is also shown in the same graph. 
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a) Sharp structures in the density of states lead to large S in simple theorites lfT9ll . 
The modern theory of the Mott transition predics a quasiparticle peak, which nar- 
rows as the transition is approached. And this could result a large thermoelectric 
response. 

b) One can think on a qualitative level of the thermoelectric coefficient as the en- 
tropy per carrier. In the incoherent regime, one could imagine that each carrier 
can transport a large amount of entropy. The incoherent regime, above a charac- 
teristic coherence temperature T*, is easy to access near a localization delcoliza- 
tion transition, because the proximity tot this boundary makes T* low. 

c) Orbital degeneracy increases the number of carriers and would be expected to in- 
crease the figure of merit. There are many orbitally degenerate three dimensional 
correlated transition metal oxide. 

Ref. |[T6l considered a model of the prototypical doped insulator LaSrTiOs, 
which has been carefully investigated experimentally [28 1. The thermoelectric prop- 
erties of this systems had not been investigated at that time. Early DMFT studies 
accounted for the divergence of the linear term of the specific heat, and the suscep- 
tibility, as well as the existence of a quasiparticle peak in the spectra |32|. 

The Hall coefficient, however, coincides with the band theory calculations, and 
is non critical near the Mott transition fl4l . It is possible to analyze the DMFT 
transport equations in two regimes: i) T <C T*, where the electronic transport is 
controlled by band-like coherent quasiparticles, well described in momentum space, 
ii) T ^> T* when the electron is better described as a particle in real space, and the 
transport is diffusive lfT6l (see below). 

The second regime is well described by the high temperature expansion, valid 
for T > D (D is the bandwidth). Approximate numerical solutions of the DMFT 
equations Ref. |16| showed that the thermoelectric response computed in the high 
temperature regime could be matched smoothly with the low temperatures results 
valid for T < T*. 



3.3 Application to LaSrTiO^ 

An approximate numerical solution of the DMFT equations for the titanides was 
shown to interpolate smoothly between the high temperature and low temperature 
region. This is consistent with the idea that DMFT reconciles the band picture at low 
energies and low temperatures, with the particle picture at high energies and high 
temperatures. The temperature scale here is set by the coherence temperature T* . 
Taking a tight binding parametrization suitable for the titanites, the thermoelectric 
figure of merit as a function of temperature and doping is reproduced in Fig. [3] 

The behavior of the thermoelectric power near the Mott transition is shown in 
Fig. [4] Notice that at low doping, the contribution from the lower Hubbard band 
dominates and the thermoelectric power is positive while at high doping the quasi- 
particle contribution dominates and the thermoelectric power is electron-like. Mea- 
surements near the Mott transition were carried out a few years later [9|, and they 
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are qualitatively, but not quantitatively, similar to the theory. This is to be expected, 
given the various approximations that were made (the electronic structure, the lattice 
distortion, and crystal field effects ignored, the impurity solvers used were very ap- 
proximate). More recent studies of these materials including lattice distortions and 
crystal field splittings have been carried out[?] but their effects on the thermoelectric 
response has not been studies. 



3.4 Low Temperature Regime 



LaSrTiOs is described by a multi-band Hubbard model. At low temperatures, the 
Fermi liquid theory is valid. The slope of the real part of the self energy scales as 
1 — 1 /Z, where Z is the quasiparticle residue. The quadratic part of the self energy 
is related to the quasiparticle lifetime, which is small in the Fermi liquid regime. 

Under these assumptions, we can rewrite a simpler expressions for the transport 
coefficients A n of a multiband Hubbard model at low temperatures 
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/: 



dx- 



x" 4>^ v (xT + ii-E'(xT)) 



cosh 2 (x/2) 



E"{xT) 



(4) 



where <P^ V is the transport function defined by <P^ V = Lk v k v k^( (B ~ e k) an d 
is the imaginary part of the electron self-energy. 

At low temperatures, Aq and A2 are simply estimated by replacing E"((o) by its 
quadratic approximation, E" (a>) • 



^{(0 2 + n 2 T 2 ) = E {2 \co). We then obtain 
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Fig. 3 Figure of merit for different values of the lattice thermal conductivity. The electron density 
is 1 —x. 
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where fio = M — ^'(0) and 



Z 2 Nk B 

2n = 7^5? *"v(«>) 



-oo4cosh 2 (x/2)[l + (x/7r) 



2U 



are numerical constants of the order unity. 

On the other hand, this approximation neglects particle-hole asymmetry and 
gives zero thermoelectricity since E\ = 0. There are two sources of particle hole 
asymmetry. One is obtained by expanding the transport function in Eq. |4]to the first 
order, which describes the particle hole asymmetry in the electronic velocities, con- 
tained in the bare band structure of the problem. This term can be approximated by 
the LDA Seebeck coefficient divided by quasiparticle renormalization amplitude Z. 
The second contribution is the result of the particle hole asymmetry of the scattering 
rate. It involves subleading cubic terms in the self energy, which scale near the Mott 
transition as 

E"(a) =Z {2) ((0)+L {3 \(0) + --- 

Z\a>)= ( - ai( ° 3+ £ a)T2) (5) 

and a\, ai are constants of order unity (even terms in frequency are not important). 
This leads to the following expression for the thermoelectric coefficient 

Ai = [^v(Mo) E\ - <£ mv (Mo) (a x El + a 2 El)/Y Q ] (6) 

where <P'(x) =d<P(x) jdx. 
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Unfortunately it has proved to be very difficult to estimate the magnitude of the 
coefficients a\ and a%. It is important to develop intuition into when these terms are 
important and their sign. Since in many cases, LDA predicts the correct sign of the 
thermoelectric power at low temperatures, perhaps the scattering time particle hole 
asymmetry Eq. (|5]l is not dominant in the LaTi03 system. This problem deserves a 
thorough investigation. 

At low temperature, the thermoelectric coefficients is 



ks k B T 



*'(w>: 



c 2 



-a 2 El 



(7) 



which clearly scales as T/Z with Z vanishing at the Mott transition. Since the linear 
term of the specific heat y scales as 1 /Z the ratio S/(yT) in a Hubbard like model 
approaches a finite value as Z vanishes: 



S 

yr 



H O(Mo) 



<£'Oo)£2 a l El + a 2 El 



<£(Mo) 



c 



(8) 



The first part of the ratio depends only on the bare band-structure quatities and is not 
effected by strong correlations. The second part, however, is due to the asymmetry 
of the quasiparticle lifetime, and might be less universal and more material and 
correlation specific. This question deserves further study. 

For the LaSrTiC>3 system, we estimated its value numerically using LDA+DMFT 
ll23l and we include its value in the plot of Behnia et.al. [4] in Fig. [2] In Ref. (4) 
it was observed that the ratio S/yT is weakly material dependent in a large number 
of materials which they compiled. From the theoretical point of view, the weak 
dependence of the ratio of Behnia et. al. on material can be view as a validation 
of the local approximation, since the most material dependence is embodied in the 
quantity Z, which cancels in the ratio S/ (yT). This suggest that the DMFT approach 
holds great promise for the search of good thermoelectric materials. Deviations from 
universality arise from the variations of the bare density of states and from the effects 
of the cubic terms in the self energy that were not included in the analysis of ref lfT6ll . 
It would be interesting to return to this problem using modern LDA+DMFT tools. 



3.5 High temperature results 

In the high temperature region, the expansion of the solution of the DMFT equations 
led to the generalized Heikes formula [30, 31] for the Seebeck coefficient. In this 
limit, thermopower is given by S = jJ./(eT), where ji is the chemical potential. 
The exact diagonalization of the atomic problem gives a set of atomic eigenvalues 
E m and their degeneracies d m . The chemical potential is then determined from the 
partition sum 

n = ^d m e-^-^/Z (9) 

m 
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where n is the number of electrons in a correlated orbital. Hence, valence of the 
solid n can be used to predict the high temperature value of thermopower. 

For the case of n < 1, which is relevant for the titanides, the expressions for 
transport quantities take the explicit form: 



a =^ D ^fi-%^ (10) 
e N{l—n) 

Here N is the spin and orbital degeneracy, and n is the electron density, D is half of 
the bare bandwidth and 70, 72 are numerical constants of order unity. 

Notice that at high temperature the Lorentz number is given by L = (k/e) 2 (D/kT) 2 y2/Yo 
with 72- Hence the Lorentz number in a model with a fixed number of particles and 
finite bandwith goes to zero at high temperatures. Thus eventually the the electronic 
thermal conductivity becomes less than the lattice conductivity and the latter con- 
trols the figure or merit. This effect was modeled in the dashed curve of Fig. [3] where 
the effects of the lattice thermal conductivity was modeled by a constant 2.0W/mK. 
The inclusion of the lattice thermal conductivity resulted in a dramatic reduction 
of the figure of merit. We can interpret the high temperature DMFT results for the 
thermal transport using a well known equation K = ivpCyZ, where vp is the Fermi 
velocity, cy the specific heat, and I the electron mean free path. Since the specific 
heat decreases as (D/T) 2 , the mean free path has saturated to a lattice spacing, and 
the velocity of the electrons is of the order of Vf. This is consistent with the value of 
the conductivity if one uses the Einstein relation a = D c dn/dfi with dn/dfi w l/T 
and the charge diffusion constant D c = Vfl. Here the mean free path I is of the or- 
der of the lattice spacing, and the Fermi velocity v/r is approximately temperature 
independent. 



4 Towards Material Design 

4.1 Rules for good correlated thermoelectricity 

From the theoretical analysis it becomes clear why Lai_ x Sr x Ti03 is not a good 
thermoelectric material near the Mott transition. The contributions from the Hub- 
bard bands and the quasiparticle peak have opposite signs, and they compete with 
each other in the interesting temperature regime, when T is comparable to T*. This 
observation leads to empirical rules for the search for good correlated thermoelectric 
materials: 
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(1) The optimal performance (when the thermal conductivity of the lattice is taken 
into account) occurs in the crossover region T mT* . Hence one should tune T* 
to the temperature region where the thermoelectric device operates. One should 
also reduce the electronic thermal conductivity (and therefore also the electric 
conductivity) until it becomes comparable to the lattice thermal conductivity, but 
not any further. 

(2) In the crossover regime, both the quasiparticle bands and the Hubbard bands con- 
tribute to the transport. Hence one should try to optimize both high temperature 
and low temperature expressions for the figure of merit. Therfore good candi- 
dates for thermoelectricity should have quasiparticle carriers and Hubbard band 
carriers of the same sign. 

We see that LaSrTi03 does not satisfy the second rule, and hence its figure of 
merit is not large. The quasiparticle contribution to the thermopower is electron-like 
while the lower Hubbard band contribution is hole like. 

It is instructive to contrast the titanites with the cobaltate materials which have a 
larger thermoelectric response. The cobaltates have holes in the lower Hubbard band 
while the quasiparticle contribution evaluated from the LDA 1 29 1 has a positive sign, 
hence it satisfies the second rule for good thermoelectricity. 

An investigation of the density driven Mott transition in the context of a two 
band Hubbard model, with one electron per site, was carried out in Ref. [22], and 
the conclusions are very similar to those obtained in the context of the doping driven 
Mott transition. 



4.2 Emergent Mottness 

Interest in thermoelectricity near the doping driven Mott transition lead to theoret- 
ical and experimental investigations of Lai_ x Sr A Ti03 and Co02Na x for small val- 
ues of the concentration parameter x. Both theory and experiment suggest that the 
thermoelectric figure of merit is not very large in this regime. On the other hand, 
the vicinity of the band insulator end, Lai-rSrrTiO^ GD and Co02Na x (see Fig. [5] 
for the phase diagram) were shown to have promising thermoelectric performance. 
Should we conclude that Mottness is bad for thermoelectricity? Not necessarily, 
after all, clear signatures of correlation were found in more realistic modeling of 
doped band insulators once the impurity potentials of the dopant atoms were taken 
into account ll26l . The impurity potential was found to restrict the spatial regions 
available for the motion of the electricity and heat carriers. In this restricted config- 
uration space, the occupancy of the electrons is close to integer and Mott physics is 
realized. 

We have suggested that similar situation occurs in the electron gas close to the 
metal insulator transition. Here, the long range Coulomb interaction generates short 
range charge crystalline lattice order. The occupancy of these lattice sites is close 
to integer filling, suggesting that the character of the metal to insulator transition 
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is that of a Wigner Mott transition [5|. The mechanism, spatial or orbital differen- 
tiation results in a restricted low energy configuration, making Mott physics rele- 
vant. This mechanism is quite general, and operates in other materials such as the 
ruthenates [20|. It could be called emergent Mottness or super-Mottness, and con- 
tains similar physics to the orbital selective Mott transition pheonomena. Hence 
(super)Mottness might be relevant for high performance thermoelectricity after all!. 
It would be useful to reconsider the most recent advances in thermoelectric materi- 
als in this light, and investigate the local magnetic susceptibility at the impurity sites 
of the high performance thermoelectrics l6l [TTI . 



5 Outlook 

The outlook for material design in the field of thermoelectric is quite promising. 
DMFT seems to capture qualitative trends in oxides of practical interest, further- 
more we have simple qualitative ideas, which can be refined and tested with tools 
of ever increasing precision. In this context, the new thermoelectric modules to be 
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developed in conjunction with the new generation of LDA+DMFT codes, look very 
appealing. In conjunction with the renewed experimental efforts in this field, the 
future looks very promising. 
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